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INTRODUCTION 


Thick-walled  cylinders  are  often  used  to  contain  very  high  pressures. ^ 

The  elastic  stress  solution  for  such  cylinders  was  developed  by  Lame^  and  is 
well  known  (see  for  example,  Reference  2).  This  solution  shows  that  for  a 
smooth  cylinder  the  maximum  stress  occurs  in  the  tangential  direction  at  the 
inside  diameter  (ID).  Such  cylinders  subjected  to  cyclic  internal  pressuriza¬ 
tion  fail  by  crack  initiation  and  growth  from  the  ID.^  To  increase  the 
fatigue  life  of  cylinders,  compressive  residual  stresses  at  the  ID  are  often 
introduced  by  autof ret tage . ^  As  a  result  of  the  autof rettage ,  a  tensile 
residual  stress  is  induced  at  the  OD.  Sometimes  after  the  autofrettage 
process,  cylinders  are  machined  on  the  OD  producing  stress  concentrations, 
such  as  keyway s ,  notches,  holes,  or  threads.  In  these  cylinders,  the  combina¬ 
tion  of  tensile  operating  stress,  tensile  residual  stress,  and  stress  concen¬ 
tration  often  result  in  the  OD,  rather  than  the  ID,  being  the  fatigue  failure 
initiation  site.^”?  OD  initiated  failures  are  less  desirable  since  they  occur 
often,  but  not  exclusively  accompanied  by  very  little  stable  crack  growth. 
Catastrophic  failure  occurs  shortly  after  the  initiation  and  coalescence  of 
many  short  penny  shaped  cracks  to  form  a  shallow  straight-fronted  crack. 6 

OD  discontinuities  in  thick  cylinders  have  been  the  subject  of  several 
experimental and  theoretical® > 10-20  studies.  The  bulk  of  the  experimental 
work5“7  >9  has  dealt  with  the  effects  of  stress  concentrations  of  fatigue  crack 
initiation,  while  one  study  has  discussed  crack  growth  from  the  OD.  All  of 
these  studies  were  on  cylinders  which  either  had  been  autof rettaged^*"®  or  were 

References  are  listed  at  the  end  of  this  report. 
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loaded  in  such  a  way  as  to  simulate  the  presence  of  autofrettage  residual 
stresses.  The  theoretical  work  has  included  internally  pressurized  cylinders 
containing  OD  cracks > 10-12 , 14 , 15 , 20  autof rettaged  cylinders  with  OD 
cracks , 16 , 19 , 20  and  autof rettaged  cylinders  with  stress  concentrations  at  the 
OD. 13, 17, 18 

In  the  experimental  studies,  Brown 5  and  Brown  and  McAlonie6>7  have 
fatigue  tested  actual  autof rettaged  cannon  barrels  containing  OD  notches. 

Their  results  showed  that  these  cylinders  failed  by  crack  initiation  at  the 
notches  followed  by  very  little  crack  growth  before  gross  through  thickness 
fracture.  In  an  attempt  to  model  this  behavior  experimentally,  without  the 
expense  of  testing  large  cannon  barrels,  Kapp  and  Underwood^  have  used  a  simu¬ 
lation  specimen.  This  specimen  was  a  modified  C-shaped  fracture  toughness 
specimen, 31  with  the  same  OD  notch  studied  by  Brown  and  McAlonie.6  By  loading 
this  specimen  properly,  the  same  stresses  were  developed  at  the  OD  notch  as 
occurred  in  the  actual  cylinder.  The  agreement  with  the  fatigue  data  in 
Reference  6  was  excellent.  By  using  this  specimen,  design  variables  such  as 
the  degree  of  autofrettage  or  notch  geometry  could  be  tested  very  efficiently. 
These  variables  were  tested  and  showed  that  fatigue  life  can  be  changed  much 
more  by  modifying  the  stress  concentration  of  the  OD  configuration,  than  by 
reducing  the  amount  of  autofrettage.  In  the  other  experimental  study,  Kapp 
and  Eisenstadt^  measured  fatigue  crack  growth  rates  in  autof rettaged  rings 
with  radial  cracks  emanating  from  the  OD.  The  results  showed  that  auto¬ 
frettage  could  increase  the  fatigue  crack  growth  rate  by  as  much  as  an  order 
of  magnitude  when  compared  to  the  results  obtained  in  specimens  which 
contained  no  autofrettage  residual  stresses. 
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The  first  theoretical  study,  by  Emery  and  Segedin,10  was  a  stress 
intensity  factor,  K,  solution  for  thick-  and  thin-walled  cylinders  loaded  by 
internal  pressure  using  a  finite  difference  method.  When  their  solutions  were 
compared  to  K  solutions  developed  using  finite  elements®  *  1®  >  ^0  and  mapping 
collocation  solutions 1^ > 16 > 19  the  agreement  was  not  good,  although  the  finite 
element  and  collocation  results  normally  agree  within  two  percent.  The 
initial  collocation  results^  were  for  cylinders  containing  from  one  to  four 
cracks  loaded  by  internal  pressure,  with  various  diameter  ratios,  W,  (W  » 
OD/ID).  The  first  finite  element  study®  was  for  a  single  crack  in  a  cylinder 
with  W  *  2. 

Subsequent  numerical  K  solutions 16 > 19 >20  included  aulofrettage  residual 
stresses  and  internal  pressure  for  as  many  as  40  standard  cracks  in  cylinders 
of  various  diameter  ratios.  Solutions  for  the  autofrettage  cases  required  the 
development  of  some  specialized  stress  analysis  techniques.  Since  auto¬ 
frettage  residual  stresses  exist  in  a  cylinder  under  the  action  of  no  external 
active  loads,  the  special  stress  analysis  methods  were  needed.  Two  basic 
methods  were  developed,  superposition  and  simulation.  Using  the  superposition 
method, >19  a  pressure  was  applied  to  the  crack  surfaces  which  was  equivalent 
to  the  negative  of  the  stress  which'  would  have  acted  at  the  same  location  in  a 
cylinder  with  no  crack  present.  In  a  modification  of  this,  Pu^  used  a  weight 
function  approach  which  was  somewhat  different  than  the  method  used  by 
Parker, 16, 19  yielding  the  same  results.  Using  the  simulation 

technique^  >19 ,20  the  autofrettage  residual  stresses  were  simulated  through  the 
use  of  an  active  thermal  distribution.  The  method,  developed  independently  by 
Parker  and  Farrow22  and  Hussain  et  al,23  showed  that  a  properly  defined,  loga- 
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rithmically  varying  temperature  distribution  produces  exactly  the  same  stress 
distribution  as  autofrettage  in  an  uncracked  cylinder.  By  applying  this  same 
temperature  distribution  to  a  cracked  cylinder,  the  K  solution  developed 
agreed  very  well  with  the  superposition  solutions.  Autofrettage  resulted  in 
very  high  K  values,  with  the  highest  values  resulting  from  the  higher  degrees 
of  autofrettage.  Also,  the  most  severe  case  occurred  when  there  were  two 
diametrically  opposed  cracks. 

The  theoretical  studies  on  notched  cylinders ^ ^ are  much  less  exten¬ 
sive  than  on  cracked  cylinders.  In  the  first  paper, ^  stresses  in  OD  notches 
were  estimated,  while  In  the  more  involved  studies , ^  the  stresses  were 
calculated  using  the  finite  element  method  and  the  thermal  loading  simulation 
method. 22, 23  The  solutions  showed  that  the  OD  notches  studied  concentrated 
the  stress  to  the  extent  that  the  notches  became  the  highest  stressed  point  in 
the  cylinder,  even  if  the  cylinder  was  not  autof rettaged .  The  autofrettage 
results  indicated  that  stresses  exceeding  the  yield  strength  of  the  cylinder 
material  acted  at  the  OD  notches.  Furthermore,  the  presence  of  the  notches 
relieved  the  autofrettage  residual  stresses  throughout  the  cylinder  by  as  much 
as  30  percent. 

In  this  report  we  will  attempt* to  relate  the  theoretical  and  experimental 
work  that  has  been  performed  to  develop  a  guideline  that  may  be  used  in  the 
design  of  thick-walled  cylinders  containing  OD  discontinuities.  The 
discussion  will  involve  three  main  areas:  special  stress  analysis  methods, 
failure  by  fatigue  crack  initiation,  and  failure  by  fatigue  crack  growth. 
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STRESS  ANALYSIS 


A  thick-walled  cylinder  of  radius  ratio  (b/a)  of  two  is  shown  in  Figure 
1 ,  which  also  defines  the  geometric  parameters  used  in  the  equations  below. 
Three  loading  conditions  will  be  discussed:  internal  pressure,  autof rettage , 
and  diametrical  line  loading  used  in  the  fatigue  crack  growth  measurements  of 
Reference  8.  The  three  loading  conditions  are  shown  schematically  in  Figure 


2. 


For  a  thick-walled  cylinder  subjected  to  internal  pressure,  the  stresses 


follow  the  well  known  Lame^  distributions 2 


(1) 


(2) 


where  p  is  the  internal  pressure  and  the  subscripts  0  and  r  signify  the 
tangential  and  radial  directions  respectively. 

For  the  case  of  autof rettage ,  the  stress  distributions  in  a  smooth 
cylinder  are  somewhat  more  complex  than  equations  (1)  and  (2).  These  auto- 
frettage  stresses  result  from  loading  the  cylinder  with  sufficient  pressure  to 
cause  plastic  deformation  to  some  extent  through  the  thickness  of  the 
cylinder.  The  degree  of  autof rettage  is  measured  by  the  elastic-plastic 
radius,  p.  This  is  the  radius  at  which,  under  the  loading  of  the  autofrettage 
pressure,  yielding  of  the  material  is  just  occurring.  The  value  of  p  gives  a 
measure  of  the  amount  of  autofrettage  defined  as  percent  overstrain  (%  over¬ 
strain  =  100  x(p-a)/t).  The  general  solutions  for  the  autofrettage  residual 
stresses  calculated  assuming  elastic,  perfectly  plastic  material  behavior  and 
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the  Tresca  yield  criteria  are:^ 

,2  n2_v2 
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(3) 


(4) 


where  o0  is  the  tensile  yield  strength.  It  is  important  to  note  that  the 
autofrettage  residual  stresses  act  in  the  cylinder  under  no  active  external 
loads.  This  causes  some  problems  when  attempting  to  determine  the  effects  of 
structure  discontinuities  on  autof rettaged  cylinders. 

Under  the  diametrical  line  loading,  the  stresses  can  be  calculated 
exactly  only  in  the  form  of  an  infinite  series. 24  The  solutions  for  thick 
cylinders  of  radius  ratio,  W,  (W  =  b/a)  ranging  from  2  to  10  have  been 
developed  by  Rigsperger  and  Davis. ^5  For  the  purpose  of  this  report,  the 
stress  which  acts  at  the  OD  for  a  cylinder  with  W  =  2  is  given.  This  stress 
acting  90°  removed  from  the  load  is: 


5.88F 

00  - -  (5) 

7Tb 


where  F  is  the  total  force  which  acts  at  the  OD  as  shown  in  Figure  2. 

To  determine  the  effects  of  OD  discontinuities  on  the  stresses  given 
above,  accurate  numerical  techniques  must  be  used  which  incorporate  such 
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discontinuities.  Three  such  Methods  have  been  used  to  do  this,  a  finite 
difference  method,  the  finite  element  method,  and  the  boundary  collocation 
mapping  technique.  The  finite  difference  approach10  was  used  to  determine 
stress  intensity  factor,  K,  solutions  for  thick-walled  cylinders  of  various  W 
with  a  single  OD  crack.  The  solutions,  when  compared  with  K  solutions 
developed  using  finite  elements8*20  and  the  collocation  method14*16*19  were 
inferior.  More  accurate  K  solutions  were  developed  using  the  latter  two 
methods.  Accurate  K  solutions  are  necessary  when  predicting  the  rate  of 
growth  of  stable  fatigue  cracks  and  unstable  fracture. 

To  develop  K  solutions  using  the  finite  element  method,  it  is  the  stress 
singularity  at  the  crack  tip  which  causes  difficulties.  A  very  large  number 
of  elements  is  usually  required  near  the  crack  tip  to  represent  the  steep 
strain  gradient.  The  recent  trend  in  finite  element  method  is  to  use  high 
order  isoparametric  elements26  so  that  a  high  degree  of  accuracy  may  be 
achieved  with  a  small  number  of  elements.  This  trend  calls  for  special 
elements27*28  which  can  produce  the  proper  order  of  singular  stresses  and 
strains  at  a  crack  tip.  The  enriched  elements27  that  contain  Ki  and  Kji  as 
additional  unknowns  are  used  to  model  the  crack  tip  region.  The  approximate 
values  of  Kj  and  Kn  are  directly  obtained  together  with  nodal  displacements. 
The  collapsed  elements,28*29  on  the  other  hand,  compute  the  nodal  displace¬ 
ments  first.  Since  a  correct  order  of  strain  singularity  is  built  in  the 
special  crack  tip  elements,  the  displacements  obtained  are  quite  accurate. 
The  stress  intensity  factors  are  indirectly  calculated  from  the  near  field 
elasticity  solution  of  displacements  in  terms  of  stress  intensity  factors. 
Both  enriched  elements  and  collapsed  elements  yield  comparably  accurate  K 


7 


values • 


In  the  modified  boundary  collocation^  method,  the  geometry  of  the 
physical  structure  is  mapped  through  the  use  of  a  complex  mapping  function  to 
a  simpler  geometry  in  a  complex  plane.  The  stress  analysis  is  then  performed 
in  the  complex  plane  and  the  solution  to  the  physical  problem  is  obtained  by 
using  the  same  complex  mapping  function  to  return  to  physical  space.  To  use 
this  method  on  cylinder  problems,  the  cracked  cylinder  is  first  mapped  to  a 
square  or  rectangular  cracked  plate.  The  plate  in  turn  is  mapped  to  a 
geometry  where  the  crack  is  represented  by  a  semicircle.  Therefore,  two 
mappings  must  be  used  to  obtain  the  solution.  This  method  was  used  to 
determine  K  for  one  to  four  OD  cracks  in  cylinders  which  had  two  W  values 
loaded  with  internal  pressure, The  K  solutions  for  one  external  crack  with 
internal  pressure  loading  are  shown  in  Table  I. 

To  use  the  stress  analysis  techniques  described  above,  the  finite  element 
model  or  the  mapped  geometry  must  be  loaded  by  active  loads.  Since  the 
autofrettage  stresses  act  in  the  absence  of  such  active  external  loads,  K 
solutions  were  developed  using  the  superposition  or  the  thermal  loading 
simulation  techniques  described  above.  Superposition  is  more  often  used  with 
the  mapping  collocation  method,  while  the  simulation  method  is  more  often  used 
with  the  finite  element  methods.  The  agreement  between  the  collocation^  and 
finite  element^  solutions  was  not  as  good  as  would  have  been  expected  (about 
five  percent  difference).  This  was  attributed  to  relatively  coarse  finite 
element  grid  used  in  Reference  20.  The  results  for  a  radius  ratio  of  two  and 
a  single  external  crack  for  autofrettage  loading  are  shown  in  Table  I. 
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TABLE  I.  K  SOLUTIONS  FOR  EXTERNALLY  CRACKED  THICK-WALLED  CYLINDERS 

(See  Table  II  for  the  definition  of  0qd) 

Internal  Pressure,  1  Crack 
K 

oOD/^ 


c/ 1 

W=2 ,  FE8 

c/ 1 

W=2 ,  COL14 

W=3,  COL14 

0.167 

1.341 

0.1 

1.20 

1.23 

0.260 

1.465 

0.2 

1.37 

1.37 

0.333 

1.635 

0.3 

1.56 

1.54 

0.417 

1.883 

0.4 

1.80 

1.82 

0.600 

2.102 

0.5 

2. 10 

2.12 

0.583 

2.404 

0.6 

2.49 

2.43  ' 

0.667 

2.753 

0.7 

- 

3.11 

0.750 

3.247 

0.833 

3.458 

Autof rettage ,  1  Crack 
K 

,  W  =  2 


c/t 

50%  Overstrain 

100%  Overstrain 

0.1 

1.129 

1.013 

0.2 

1.235 

0.977 

0.3 

1.424 

0.969 

0.4 

1.647 

0.5 

1.924 

0.6 

2.041 

0.923 

9 


To  use  the  numerical  solutions  given  in  Table  I  to  predict  the  rate  of 
growth  of  fatigue  cracks  or  fracture,  K  should  be  known  as  a  continuous 
function  of  crack  length.  This  can  be  accomplished  by  fitting  a  polynomial  to 
the  numerical  K  solutions  along  with  limiting  solutions. 30  The  expressions 
that  were  developed  for  some  of  the  numerical  results  in  Table  I  have  the  form 

K  =  oqd  »^c(1.12  +  A^(c/t)  +  A2(c/t)z  +  +  A4(c/t)4) 

where  oqd  is  the  stress  which  acts  at  the  OD  of  an  uncracked  cylinder  due  to 
the  applied  loading;  ,  A2 ,  A3  and  A4  are  constants;  c  is  crack  length;  and  t 
is  wall  thickness  (b-a).  The  factor  1.12  arises  from  the  fact  that  as  c 
approaches  zero,  K  approaches  the  following :^1 

K  =  1.12  oqd  /ttc  ,  (6) 

The  values  of  the  constants  for  various  loadings  are  given  in  Table  II 
for  the  cases  of  internal  pressure,  diametrical  line  loading,  50  percent 
overstrain  and  100  percent  overstrain,  for  a  cylinder  with  W  =  2  and  a  single 
external  crack.  These  expressions  represent  the  numerical  results  within 
three  percent  for  any  of  the  loadings  considered. 

As  stated  previously,  there  have  been  fewer  theoretical  studies  on  thick 
cylinders  which  contain  OD  discontinuities,  such  as  notches,  rather  than 
cracks  on  the  OD.  In  these  studies,  the  stresses  have  been  calculated  by 
either  estimation^  or  with  finite  elements.  12,18  The  most  extensive  work^ 
studied  the  case  of  a  thick  cylinder  with  W  =  1.74,  subjected  to  either 
internal  pressure,  100  percent  overstrain  or  60  percent  overstrain.  The 
cylinder  contained  a  single  0D  notch  which  had  a  root  radius  to  thickness 
ratio  (R/t)  of  0.013,  and  a  depth  to  thickness  ratio  (d/t)  of  either  0.088  or 
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TABLE  II.  COEFFICIENTS  FOR  THE  EXPRESSIONS  FOR  THE  STRESS  INTENSITY  FACTORS 
K  =  O0D(^*c(1* 12  +  Ai(c/t)  +  A2(c/t)2  +  A3(c/t)3+  A4(c/t)4) 


Loading  Condition 

Al 

a2 

a3 

a4 

Internal  Pressure® 

0.31 

6.85 

-12.12 

+10.02 

Diametrical  Lineload® 

-1.14 

1.23 

-  2.20 

0.0 

50%  Overstrain 

0.0 

0.0 

16.9 

-21.0 

100%  Overstrain 

-1.44 

+4.444 

-  4.31 

0.0 

Loading  Condition 

Internal  Pressure 


Diametrical  Load 


50% 


100% 


°OD 

2pa2 

b2-a2 


5.88F 

irb 

2oo  x  lx 

/3  4  w2-l  4 


2o0 

/3 


2Jln  w 

(1  +  - - ) 

w2-l 


w+1 

2ln  ( - )) ,  x  =  (1  + 

2 


w 
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0.177.  The  results  of  the  finite  element  calculations  are  given  in  Table  ICE. 
These  show  that  the  notches  constitute  a  significant  stress  concentration 
which  produced  a  stress  concentration  factor  kt  as  high  as  6.6  in  the  case  of 
the  deep  notch,  100  percent  overstrain  condition.  Another  interesting  finding 
of  the  study  was  that  the  presence  of  the  notch  caused  a  reduction  or  a 
redistribution  of  the  autof rettaged  residual  stresses  throughout  the  cylinder. 
Therefore,  also  included  in  Table  III  are  the  relief  factors  Rg  ,  which  are  the 
ratios  of  the  remaining  residual  stress  to  the  residual  stresses  which  acted 
in  the  cylinder  prior  to  the  introduction  of  the  notch.  To  use  the  factors  in 
Table  III  to  calculate  the  stresses  which  act  at  the  OD  notch  the  following 
expression  is  used: 

anotch  =  *s  a0(r  =  6— d)  (7) 

where  o@(r)  is  given  by  equation  (3). 

TABLE  III.  STRESS  CONCENTRATION  FACTORS  kt  AND  STRESS  RELIEF  FACTORS  kg 
FOR  AN  OD  NOTCHED  CYLINDER,  W  =  1.74,  ROOT  RADIUS  OF  THE 
NOTCH  =  0.013  t 


Loading  Condition 

Deep  Notch 
(Depth-0.178  t) 

Shallow 
( Depth=0 

Notch 
.088  t) 

kt 

Rs 

kt 

*s 

Internal  Pressure 

5.4 

- 

3.5 

- 

60%  Overstrain 

5.8 

0.71 

3.5 

0.72 

100%  Overstrain 

6.6 

0.70 

3.7 

0.70 
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PREDICTING  THE  FAILURE  OF  CYLINDERS  WITH  OD  DISCONTINUITIES 


Using  the  results  from  the  various  stress  analyses  described  above,  we 
will  attempt  to  predict  failure  under  fatigue  conditions.  These  failures  can 
occur  by  either  of  two  mechanisms,  fatigue  crack  initiation,  or  fatigue  crack 
growth.  In  the  case  of  crack  initiation,  three  factors  contribute  to  failure, 
the  alternating  stress  from  external  loading,  tensile  mean  stress  from  the 
autofrettage  stresses,  and  the  stress  concentration  from  the  OD  discontinuity. 
For  failure  by  fatigue  crack  growth,  two  factors  must  be  considered,  the 
alternating  K  from  the  applied  loading  and  the  ratio  of  minimum  K  (Kmin)  to 
maximum  K  C^^),  during  fatigue  due  to  changes  in  autofrettage  conditions. 
Crack  initiation  is  considered  first. 

To  predict  fatigue  crack  initiation  a  parameter  must  be  developed  which 
includes  all  three  of  the  contributing  factors  described  above.  Fuchs ^2  has 
proposed  such  a  parameter  for  predicting  crack  initiation  at  long  life  ( > 1 0 ^ 
cycles).  He  states  that  if  the  following  equation  is  valid  crack  initiation 
occurs : 

kf  Aa 

ae  <  —  (—  +  Cam)  (8) 

n  2 

where  ae  is  the  endurance  limit  of  the  material,  kf  is  the  fatigue  stress 
concentration  factor,  Aa  is  the  range  of  stress  applied  (Aa  =  amax  -  ^min) >  C 
is  constant  with  a  value  of  about  0.5,  and  ara  is  the  mean  stress  (am  =  (crmax  + 
amin)/ 2) . 
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Extending  Fuch’s^  work  to  short  life  situations  (<  106  cycles),  it  is 
assumed  that  the  value  of  stress  predicted  by  the  right-hand  side  of  equation 
(8)  is  equivalent  to  the  uniaxial  alternating  stress  which  results  in  less 
than  10^  cycles.  The  relationship  between  the  stress  and  life  is  determined 
by  the  S-N  property  of  the  material  studied.  The  S-N  behavior  of  a  steel  can 
be  estimated  to  a  first  order  approximation  using  the  following  assumptions . 33 
At  a  life  of  103  cycles,  a  specimen  in  a  compression-tension  fatigue  test  will 
fail  with  an  alternating  stress  of 

A  0 

oa  =  —  =  0.9  outs  (@  N  =  103  cycles)  (9) 

2 

where  auts  is  the  ultimate  tensile  strength  of  the  material.  At  a  life  of  106 
cycles,  the  allowable  alternating  stress  is  approximated  as 

Aa 

°a  **  =0.5  kakbkc  ^uts  (10) 

2 


where  ka  is  a  factor  related  to  surface  finish,  k^  is  a  factor  related  to  the 
size  of  the  component  considered,  and  kc  is  a  reliability  factor  dealing  with 
statistical  considerations.  Assuming  the  steels  considered  have  an  ultimate 
tensile  strength  of  about  1250  MPa,  the  cylinders  have  a  machined  finish  and 
are  larger  than  2  cm,  and  we  desire  99  percent  reliability,  the  respective 
values  of  ka,  k^,  and  are  0.65,  0.85,  and  0.868.  Equation  (10)  then 
becomes 


a 


a  " 


Aa 

2 


0.24  auts  (@  N  =  106  cycles) 


(ID 


It  is  also  assumed  that  the  S-N  property  of  the  material  falls  on  the 
straight  line  in  log  S-log  N  space  which  connects  the  two  points  defined  by 
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equations  (9)  and  (11).  The  solid  line  in  Figure  3  is  the  estimated  S-N 
property  of  the  material.  The  data  points  in  the  same  figure  are  actual 
measurements  of  fatigue  life  from  References  8  and  9  using  Fuch’s  criterion 
(Equation  (8)).  The  stress  concentration  factors  for  the  data  from  Reference 
9  are  from  Reference  17  and  are  considered  to  be  quite  accurate.  For  the  data 
from  Reference  8,  the  stress  concentration  factors  were  estimated  using  a 
Nuber  diagram,  34  and  are  considered  less  accurate  than  those  calculated  in 
Reference  17. 

The  comparison  of  the  available  data  with  the  estimated  fatigue  life  to 
crack  initiation  suggests  that  the  procedure  outlined  above  is  conservative  if 
the  stress  concentration  factors  are  known  accurately.  If  the  stress 
concentration  factors  must  be  estimated,  the  failure  prediction  gives  an 
approximate  average  of  failure.  In  some  cases  life  is  underestimated,  in 
others  it  is  overestimated.  To  use  this  method  in  design  warrants  the  use  of 
safety  factors  if  kt  is  not  accurately  known.  The  dotted  line  in  Figure  4 
shows  the  failure  prediction  using  a  safety  factor  of  two  on  stress.  In  this 
instance,  all  of  the  data  fall  well  beyond  the  predicted  failure  curve  when  kj- 
is  estimated.  Using  safety  factors  should  only  be  applied  when  kt  is 
estimated,  since  they  may  be  overly  ‘restrictive  when  kt  is  accurately  known. 

To  predict  failure  by  fatigue  crack  growth,  a  crack  growth  law  which  is  a 
property  of  the  material  is  needed*  The  first  law  developed  was  the  well 
known  Paris  power  law^5 

dc 

—  “  Cj  AK™  (12) 

dN 
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where  C\  and  m  are  constants,  dc/dN  is  the  measured  fatigue  crack  growth  rate, 
and  AK  is  the  range  of  stress  intensity  factor  applied  during  loading.  This 
law  is  useful  in  predicting  crack  growth  if  K  ranges  from  zero  to  a  positive 
value.  If  K  varies  from  a  positive  value,  K^n,  to  a  maximum  value,  K^x ,  the 
crack  growth  rate  may  be  different  from  that  predicted  using  equation  (12). 
This  is  the  case  when  considering  the  effect  of  autofrettage  on  crack  growth. 8 
To  model  the  effects  of  cycling  between  two  positive  K  values,  several 
laws  have  been  suggested.  These  laws  have  also  incorporated  static  properties 
of  the  material  in  their  formulation.  Either  of  two  static  properties, 
fracture  toughness,  Kxc ,  or  yield  strength,  aQ  are  used.  A  crack  growth  law 
using  Kxc  is  of  the  form^ 

dc  C2 AKm 

—  - -  (13) 

dN  {(l-Rf)Kxc  -  AK> 

where  C2  and  m  are  constants  and  Rf  is  the  ratio  of  minimum  K  to  maximum  K 


( ^f  =  ^min/^max)* 

An  example  of  a  crack  growth  law  incorporating  the  yield  strength  of  the 
material  is  due  to  Tabeshfar  and  Williams :  37 


dc 

—  —  sa 

dN  . 


v  tn 
^ax 

C3( - ) 


(14) 


where  C3  and  m  are  constants. 

We  used  these  laws  to  predict  crack  growth  through  the  autof rettaged 
cylinders  reported  in  Reference  8.  The  material  used  in  that  study  was  an 
AISI  4340  steel  with  yield  strength  of  about  1140  MPA,  and  Kxc  of  about  150 
MPa/m.  The  fatigue  crack  growth  property  (Equation  (12))  has  been  measured 
for  this  material. 38  Using  these  data  the  constants  in  all  the  crack  growth 
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laws  can  be  determined.  For  the  Paris,35  Forman  et  al,36  and  the  Tabeshfar 


and  Williams3^  models  respectively,  the  crack  growth  laws  are: 

dc 

—  =  2.46  x  10-b  AK2*! 
dN 

dc  4.80  x  10-8  AK2 
dN  { 150( 1-Rf )  -  AK} 


(15) 

(16) 


,u  2.71  x  10  Kmax^’^ 
dN 


(17) 


All  of  these  growth  laws  are  variable  separable  differential  equations  of 
the  form 


dc 

7"  =  F(K,Rf,...)  (18) 


To  make  crack  growth  predictions  the  crack  growth  laws  must  be  integrated  to 
give  a  function  which  relates  crack  length  to  number  of  cycles,  N.  This  is 
accomplished  by  rearranging  equation  (18)  to 


dN 


dc 

F( K,Rf , . . . ) 


(19) 


Equation  (19)  is  now  in  a  form  which  can  be  easily  integrated  using  numerical 
techniques,  the  simplest  being  the  trapezoidal  rule.3^ 

In  Reference  8  crack  growth  measurements  were  obtained  in  cylinders 
which  were  not  autof rettaged  or  autof rettaged  to  either  50  percent  or  100 
percent  overstrain.  All  three  crack  growth  laws  were  used  to  predict  crack 
growth  under  all  three  of  the  imposed  testing  conditions,  using  the  stress 
intensity  factor  expressions  from  Table  II  to  calculate  the  required  K 
relationships  for  the  predictions.  The  results  are  shown  in  Figures  4,  5,  and 
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6  for  the  non-autof rettaged ,  50  percent  overstrain  and  100  percent  overstrain 
conditions  respectively. 

For  the  case  of  no  autof rettage ,  Figure  4,  it  is  clear  that  all  three 
crack  growth  models  adequately  model  the  reported  crack  growth  from  Reference 
8.  There  is  an  ordering  of  the  crack  growth  predictions  which  shows  that 
the  Paris  model  predicts  the  slowest  crack  growth,  the  Tabeshfar  and  Williams 
the  fastest  and  the  Forman  et  al  falling  in  between  giving  the  best  approxima¬ 
tion  to  the  actual  behavior.  This  suggests  that  there  is  an  effect  of  Rf 
ratio  even  in  the  case  of  no  autofrettage  when  one  would  consider  the  Paris 
law  to  provide  the  best  prediction.  In  the  experiments  reported  in  Reference 
8,  Rf  was  0.1,  while  the  crack  growth  data  from  Reference  7  was  developed  with 
Rf  =  0. 

In  examining  the  results  in  autofrettage  cylinders,  Figures  5  and  6,  we 
find  that  no  model  gives  an  accurate  estimate  of  fatigue  crack  growth.  The 
poorest  correlation  occurs  in  the  case  of  100  percent  overstrain,  where  the 
Tabeshfar  and  Williams  model  predicts  the  crack  should  grow  to  a  length  of 
about  9.5  mm  in  less  than  2,500  cycles  when  it  actually  required  about  50,000 
cycles  to  obtain  that  depth.  To  explain  this  difference  it  should  be  noted 
that  in  the  initial  crack  growth  regions,  the  Paris  model  gives  a  relatively 
good  approximation  for  the  crack  growth  behavior.  At  greater  crack  depths, 
the  Paris  prediction  substantially  overestimates  the  cycles  required  to  grow 
the  crack  further.  Although  the  Forman  et  al,  and  Tabeshfar  and  Williams 
models  underestimate  the  crack  growth  behavior  initially,  the  cycles  required 
to  cause  additional  crack  growth  are  reasonably  approximated.  For  example 
consider  the  50  percent  overstrain  case,  Figure  5.  To  grow  a  crack  from  4  mm 
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to  about  7.5  mm  requires  about  15,000  cycles  according  to  Forman  et  al.  In 
the  data  from  Reference  8  for  the  same  range  of  crack  growth,  25,000  to  30,000 
cycles  were  actually  needed.  These  observations  suggest  that  near  the  starter 
notch  used  to  initiate  the  fatigue  cracks,  a  model  assuming  no  autofrettage 
applies ,  while  removed  from  the  starter  notch,  a  model  which  incorporates  a 
significant  contribution  from  autofrettage  seems  to  give  a  reasonable  estimate 
This  indicates  that  the  presence  of  the  notch  may  significantly 
change  the  residual  stresses  near  the  notch.  Removed  from  the  notch,  there 
seems  to  be  little  effect  on  the  residual  stress  due  to  the  notch;  the  crack 
growth  model  of  Forman  et  al  gives  a  reasonable  estimate  of  crack  growth. 

CONCLUSIONS 

Much  work,  mostly  theoretical  stress  analyses,  has  been  performed  to 
study  the  effects  of  OD  initiated  failures  in  thick-walled  cylinders.  The 
stress  analyses  were  developed  using  very  accurate  numerical  techniques  which 
have  the  potential  for  being  very  useful  in  the  design  of  thick-walled 
pressure  vessels.  At  this  time  there  seems  to  be  insufficient  experimental 
results  to  use  the  theoretical  studies  to  predict  the  OD  initiated  failure  of 
these  vessels  with  any  great  accuracy.  Using  the  analytical  results  and 
estimated  fatigue  crack  initiation  properties,  conservative  estimates  of 
atigue  life  can  be  made  which  may  be  useful  to  prepare  cylinder  designs  which 
should  not  fail  by  fatigue  crack  initiation.  Also,  if  failure  occurs  by 
fatigue  crack  growth,  conservative  estimates  of  fatigue  life  may  be  obtained 
using  available  crack  growth  models.  Neither  of  the  two  models  which  account 
or  the  mean  stress  affect  on  fatigue  crack  growth  agreed  very  well  with  the 
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observed  data.  The  use  of  the  model  by  Forman  et  al^6  did  give  the  best 
estimate  of  the  maximum  fatigue  crack  growth  rates  and  should  give  a 
conservative  estimate  of  fatigue  life.  More  work  is  needed  to  develop  the 
technology  necessary  to  estimate  fatigue  failure  with  greater  accuracy. 
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Figure  1. 


A  Thick-Walled  Cylinder  Containing  Various  OD  Discontinuities. 
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Figure  2.  The  Three  Loading  Conditions  Considered.  . 
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Figure  3.  Comparison  of  Measured  Fatigue  Life  and  the  Estimated  Fatigue 
Life  Expanding  on  Fuch's  Crack  Initiation  Criteria. 32 
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Figure  4.  Comparison  of  Measured  Fatigue  Crack  Growth  and  Predicted 
Fatigue  Crack  Growth  Using  Various  Models  for  a  Cylinder 
Containing  No  Residual  Stress. 
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Figure  5.  Comparison  of  Measured  Fatigue  Crack  Growth  and  Predicted 
Fatigue  Crack  Growth  Using  Various  Models  for  a  Cylinder 
Containing  50  Percent  Overstrain  Residual  Stresses. 
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Figure  6.  Comparison  of  Measured  Fatigue  Crack  Growth  and  Predicted 
Fatigue  Crack  Growth  Using  Various  Models  for  a  Cylinder 
Containing  100  Percent  Overstrain  Residual  Stress. 
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